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The magic values for B m and D. m rely on the precise determination of the atomic parameters for the optical clock 
transition under investigation. In order to derive the values reported within the text of the present work, we have 
discovered that the published information on the required atomic parameters is some time confusing. The main 
target of this Supplemental Material is to make clear some notations, and to list the numerical values of the Zeeman 
parameters used in the calculation in order to provide a correct comparison between the dressing of the different 
optical clocks. In addition we precise our definition of the rf parameters in order to derive the precise value of the rf 
amplitude to be applied in the optical clock experiments. 

We define the Rabi angular frequency Q,j as: 

0, = (1) 

where B r f is the amplitude of the oscillating rf field, gj the Lande factor of the j = (S.P) level, fip the Bohr 
magneton. This angular frequency is equivalent to the Larmor frequency Wi defined in [l| and associated to the 
classical rf field amplitude. Our definition is four times larger than the frequency b introduced in 0, S| • This explains 
why the coefficients of the function S(x) of Eq. (3) in the main text arc different of those of Q. 

The linear nuclear spin dependent Zeeman shift on 3 Pq and 1 Sq states depend on the Lande factors gs and gp, 
see Eq. (1) of the main text. Their values used in the numerical calculations are listed in Table [H For the alkaline- 
earth-metal atoms, 87 Sr and 171 Yb, gs is calculated from the value of the shielded nuclear magnetic moment Jj]. The 
gp values, gp = gs + Sg, are obtained from the experimental values of 5g, for 87 Sr in Q and for 171 Yb in @ . The 
coefficient Av B of the quadratic term was experimentaly measured for 87 Sr in 0]. That term has been measured on 
the fermion 171 Yb [f|, but we use the more precise boson 174 Yb value of 0], since the quadratic Zeeman shifts of both 
isotopes should be identical, as measured for the boson M an d fermion pj Sr isotopes. For 199 Hg the gs factor is 



calculated from the value of the shielded nuclear moment [10. Illj. gp was measured in 12|, A.v B is given in [1 



TABLE I: Atomic parameters of 1 So and 3 Po fermionic states used in the main text 



Atom 


87 Sr 


171 Yb 


199 Hg 


I 


9/2 


1/2 


1/2 


gs^B/h (Hz/juT ) 


1.844 [4J 


-7.439 [4] 


-7.590 [10] 


S gl i B /h (Hz/jiiT) 


1.084(4) [5J 


-4.20[6J 


-5.712 


g P ^ B /h (Rz/fiT) 


2.928(4) 


-11.64 


-13.30 [12] 


gp/gs ratio 


1.59 


1.56 


1.75 


A^/h (Hz/mT 2 ) 


-23.5(2) [7] 


-6.1(1) [8] 


-2.44 [13] 



Note that the numerical values of Tabic I of the main text critically depend on the numerical values used for gs and 
gp (not Sg only), and for . The comparison between the strontium, ytterbium and mercury parameters evidences 
that the ratio of the nuclear magnetons between excited and ground states is similar for the three species. Instead 
a large difference is associated to the second order Zeeman shift, being larger in strontium by a factor between four 
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and ten. Such difference leads to the magic magnetic field values reported in the Table I of the main text at fixed oj r f 
very low for ytterbium and mercury. For an implementation of the dressing of those atoms, it may be convenient to 
increase u> r f and bring the magnetic fields in a range where a better control is reached. 
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A Zeeman-insensitive optical clock atomic transition is engineered when nuclear spins are dressed 
by a non resonant radio-frequency field. For fermionic species as 87 Sr, 171 Yb, and 199 Hg, particular 
ratios between the radiofrequency driving amplitude and frequency lead to "magic" magnetic values 
where a net cancelation of the Zeeman clock shift and a complete reduction of first order magnetic 
variations are produced within a relative uncertainty below the 10 _1S level. An Autler-Townes 
continued fraction describing a semi-classical radio-frequency dressed spin is numerically computed 
and compared to an analytical quantum description including higher order magnetic field corrections 
to the dressed energies. 

PACS numbers: 32.80. Qk, 06.30. Ft, 32.30.Bv 



The interaction between magnetic fields and the nu- 
clear/electronic magnetic moments represents a flexible 
tool for the control of the internal and external degrees of 
freedom in atoms or molecules, widely employed in pre- 
cision measurements, frequency metrology, and coherent 
manipulations of quantum systems. In frequency metrol- 
ogy, the presence of magnetic fields may represent a limit 
on the realization of specific targets. In the attractive 
context of the " magic" wavelength combining a vanishing 
differential shift of the clock levels with the Lamb-Dicke 
regime greatly reducing the motional effects [H-IH, the 
quest for "magic" magnetic field values where first or- 
der Zeeman shift and magnetic fluctuation of the atomic 
transition is annulled, was proposed in Q and studied 
experimentally in rubidium [Bj]. Best performances in 
optical clocks are accessible by using atomic transitions 
allowed by a weak hyperfinc mixing mediated through a 
small spin-orbit coupling with a resolution at the mHz 
level. To obtain even better performances, it has been 
proposed to use the bosonic even isotopes eliminating 
the nuclear spin and removing completely the first order 
Zeeman effect with a residual second order magnetic shift 
comparable to those of the ion standards . However in 
order to get rid of strong cold collision frequency shifts as- 
sociated to the bosons, the present frequency metro logy 
is concentrated on the fermionic species, as 87 Sr |7H10l|. 
171 Yb [HEl, 199H S [3 El- There, the first and second 
order Zeeman shifts contribute by one order of magnitude 
above the projected 10~ 18 fractional uncertainty of the 
frequency standard @, H, • 

This work focuses on the atomic magnetic moment en- 
gineering with the target of getting rid of atomic proper- 
ties sensitive to external electromagnetic fields. A scheme 
based on the radiofrequency (rf) quantum engineering of 
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fermionic atomic states is presented in order to produce 
levels experiencing a vanishing first-order Zeeman clock 
shift. The cancelation of first-order Zeeman shift applies 
also to the vectorial ac Stark shift, equivalent to an ef- 
fective magnetic field, the only contribution of this kind 
appearing in 171 Yb and 199 Hg. The basic idea of let- 
ting a paramagnetic system mimic a non-magnetic one 
originates from artificial or synthetic magnetism where 
an atomic Hamiltonian is created by proper electromag- 
netic fields in order to simulate a given magnetic config- 
uration [l5j . Our work is inspired by the dressed-atom 
rf quantum engineering fl6| , where the paramagnetic re- 
sponse for two species, atoms in [l~7j or atom/neutron 
in [Tc| . is tuned into the resonance. 

The cancelation of the first order Zeeman effect is pro- 
duced by the atomic dressing at a rf frequency much 
larger than the effective Larmor precession, equivalent 
to a frequency modulation of the nuclear magnetization 
and a shielded nuclear response to the static magnetic 
field. The different rf response for the ground and ex- 
cited states of the clock transition leads to crossing nodes 
in the energy diagram, where the atoms become non- 
magnetic. This change from a paramagnetic system to an 
non-magnetic one shares a strong analogy with a Landau 
theory of phase transition. In addition at magic static 
field values the rf dressing engineers a Zeeman-insensitive 
atomic clock. This magic cancelation arises from the non- 
linear magnetic Hamiltonian associated to the rf dressing 
of the two-electron system. Even if the dressing does not 
eliminate the second-order Zeeman contribution, its con- 
tribution to clock state separation is strongly decreased 
by an operation at a magic magnetic field. The stability 
of the ratio between rf dressing amplitude and rf angu- 
lar frequency required to produce a target non-magnetic 
state matched to the aimed optical clock accuracy is ex- 
perimentally reachable. The present approach of a magic 
magnetic field cannot be applied to tensorial ac shifts. 

For the alkali-earth atoms, a modified Breit- Wills the- 
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FIG. 1: (Color online) (a) Bare Zeeman energy splitting E mF 
of the 87 Sr ^S^rnp = ±9/2) and | 3 P ;m F = ±9/2) clock 
states versus magnetic field B. (b) Dressed Zeeman energies 
£™ F ( 3 Po, 3 S ) vs x = ftyw r/ at o;/2tt = 2 kHz and B=25 
fiT. The crossing nodes with a zero first order Zeeman shift 
determine the magic rf values. 



energies of the highest |?tif| 87 Sr clock levels vs B are 
plotted in Fig. 1(a). 

We drive the clock atoms by a non- resonant rf field, lin- 
early polarized and orthogonal to B, at angular frequency 
w r f and (ground-state) Rabi frequency Q, see Supplemen- 
tal Material. A strong modification of the Lande g-factor 
occurs in the regime where (J-bB <C hu r f. For dressing 
by a large number of rf photons, a pcrturbativc quan- 
tum analysis predicts a dressed Lande g-factor depen- 
dent on the zeroth-order Bessel function of the first kind 
9j(xj) = g±Ip(xj), with Xj = Q.j/u] r f, % = gjtt/gs, 
(j = P, S) [21| . That dependence was verified in experi- 
ments on atoms nam 1221. n eutrons [23| and chromium 
Bose-Einstein condensate [24| . It is valid for whatever 
spin value and equally spaced Zeeman levels f25j . When 
HbB ~ hu! r f, the g'j expression includes an additional B 
dependence given by [25|, [26[ 



d 
9j 



(Jj 



J (Xj) 



( 9jHBB \' 



S(xj 



(2) 



S(x) being a product of Bessel functions. However for the 
first two crossing nodes of the rf dressed Zeeman energies 
of Fig. 1(b) the following approximated analytical S(x) 
expression given by [2tJ provides the required accuracy: 



S(x) 



16 



2025a; 4 



[a(x)J 2 (x) + P(x)Ja(x) - 7(2) Jq(x)] 



(3) 

where functions are a(x) = 75 (5a; 2 — x 4 /4) , (3(x) = 
6 (408 - 74a; 2 - 23x 4 /16) , y{x) = 145 (3a; 2 - a; 4 /2) /49. 
When the dressed Lande g-factor of Eq. @ is substituted 
into the energies of Eq. ([I]) , the rf dressed energies con- 
tain both B 2 and B 3 nonlinear terms. 

We derive the exact rf dressed Zeeman energies 
Ef np (j) from the Autler-Townes continued fraction 



ory describes the action of a magnetic field B produc- 
ing linear and quadratic nuclear spin dependent Zeeman 
shifts for the doubly forbidden | 1 So) — > | 3 Po) optical 
clock transition [h|. For the \F,mp > Zeeman level 
the energy E mp is 



E P ,m p , ( 3 Po) = m Fl gpn B B + A^B 2 , 
Es,m F CS ) = m F gsiJ-BB, 



(1) 



where m' F and m F are the lower/upper magnetic quan- 
tum numbers, gp and gs the Lande g-factors, with 

(2) 

9p = gs + ^5; and the second order Zeeman con- 
tribution, mp independent. We will focus our attention 
on fermionic systems spin polarized in the extreme Zee- 
man sub-levels fl9l . |20| where a systematic average on 
the transition frequencies of optical transitions symmet- 
rically placed around line center is currently applied to 
cancel the linear Zeeman shift and to probe accurately 
the second order Zeeman correction. Their parameters 
are listed into the Supplemental Material. The Zeeman 



E d mF {j) = E mF (j) + 2m F MlZL • L(j). 

gs 



(4) 



The function L, representing the m = 1/2 eigenvalue, 
normalized to the Rabi frequency [2l|, for a spin- 1/2 sys- 
tem having W21 energy splitting and rf dressed by a u> r f 
field with Rabi frequency {I, is given by [28[ 
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"21 (i_ OM 



L-2- 
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The dressed nuclear Lande g-factor is obtained deriving 
the dressed energies with respect to B. 

The dressing field strongly modifies the Zeeman ener- 
gies for each thf clock state as shown in Fig. QJb) for a 
given static field B. The energies follow mainly the zero- 
order Bessel function dependence. The different l S$ and 
3 Po sensitivity to the dressing created by Sg produces 
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FIG. 2: (Color online) Dressed Zeeman shifts of the 87 Sr 
rriF = —9/2 or ulf = 9/2, in alternance, n clock transi- 
tion vs B for u r f /27T = 2 kHz and the fi/cj r / values at, or 
around, the crossing nodes in Fig. [T] The curvature around 
the B values of the shift minima is the dressed second order 
Zeeman shift. For the (b) open red data the x is modified 
by one part in 10 6 from the crossing value to a value where 
the second order Zeeman shift is annulled, in presence of con- 
stant Zeeman bias. Open and closed dots • based on Eqs. Q 
and (JSJ. Lines based on Eqs. and Q are approximated 
solutions providing a good description around the first two 
crossing nodes only. 10 /iHz correspond to a 2.10 -20 clock 
fractional shift. 

several crossings between the clock energies at particular 
Q,/uj r f ratios. The dressed Zeeman clock shift 

A< F ^„ 4 - Ef n , F ( 3 P ) - E d mp ('So) (6) 

is exactly compensated for specific rf dressing param- 
eters, for instance in Fig. (5Ja) at the x = il/u> r f w 
0.96, 3.11, 5.44, ... for a tt transition. The oj r f and fl com- 
pensating values are determined by imposing the clock 
transition to be immune from the Zeeman shift, i.e., equal 
dressed magnetic energies. 

A more ambitious target is to derive the " magic" fields 
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FIG. 3: (Color online) (B,x) branches, derived from Eqs. ©, 
defining a zero clock Zeeman shift AE d = (solid red line) 
and a zero derivative dAE d /dB = (dashed-dot red line) for 
the alternating vtlf = —9/2 and tuf = 9/2 n clock transitions 
in 87 Sr. Their intersections determine the (B m ,x m ) magic 
values also denoted by the solid dots • from a numerical eval- 
uation of Eq. and Eq. ©. The continuous line dependence 
reports the B(x) of Eq. ([8]), but magic B m values occur only 
at magic x m parameters. All curves for iu r f/2n = 2 kHz. 

where the pure-Zeeman differential shift is zero and also 
independent of the field value. That applies to the 
dressed Zeeman clock shift plotted in Fig. [Ha) with os- 
cillations in the B dependence. At w r //27r = 2 kHz and 
Q/w r f = 0.9638744 we have a magic magnetic field value 
B m w 2.06 /iT which exactly cancels the full Zeeman 
shift for the mp = —9/2, tt transition, experiencing in 
addition a reduced (B — B m ) 2 sensitivity for the Zeeman 
energies. A reduction by one-order of magnitude for that 
second order sensitivity is obtained operating close to 
the higher order crossing nodes of Fig. [IJb), the first two 
also shown in Fig. [Ha) and corresponding to alternating 
mp = —9/2 and mp = 9/2 Zeeman sub-levels. Fig. [^b) 
shows that, by changing the x parameter, a perfect can- 
celation of the (B — B m ) 2 magnetic sensitivity is reached, 
at expense of a clock constant bias. 

The "magic" (B m ,x m ) values where the Zeeman shift 
and the first order sensitivity to weak field variations are 
simultaneously canceled are derived by imposing 

AE d mF _ mF ,(x)=0; ^ Hr^J = ' ( ? ) 

The continued fraction solution of Eqs. (|4]) and ([5]) de- 
termines the magic values associated to the above con- 
ditions. Fig. [3] reports the graphical approach applied 
to derive these magic B m values. Notice that a suppres- 
sion cannot be realized simultaneously for more than one 
mp spin dependent transition because the B square de- 
pendence of Eq. @ imposes for each optical transition a 
matched dressed g-factor compensation. 



4 



TABLE I: Ma : 
199 Hg optical 
fixed uj r f /2ty - 
with the 10" 6 
below 10" 19 . 



_ic (x m , B m ) pairs for it and a S7 Sr, 171 Yb and 
clock transitions based on Eqs. ((4} and (O, at 
= 2 kHz. The x m = 0. m /uj r f values are reported 
accuracy required to produce a fractional shift 







87 Sr 




rriF — > vfipi 


_9 _j, _9 


^2 ^ ^2 


— - — > — - +- — > + - 




0.963875 


3.110838 


5.444392 7.776908 


B m (fiT) 


2.06 


5.1 

171 -i^i 


9.0 20.0 



rriF rripi 

B m OT) 
Tn p — ¥ rapi 

B m (flT) 



+ 1 -> +\ 

0.977571 

0.08 
-i -> +i 
1.825873 



—I -> -i 
3.157095 
0.12 



+ 2" — *■ ~~ 

4.107330 



+ 1 -> +1 —I — > —I 
5.527470 7.905770 
0.33 0.68 

5.542963 6.953892 2 



tuf — ► "If' 

Bm (JUT) 

TTiF —5- mpi 



0.11 


0.59 


0.88 0.82 




lyy Hg 




+1 -> +1 

0.911529 


—I -> —I 


+1 +1 —I — > —I 


2.930862 


5.116966 7.220859 2 


0.02 


0.05 


0.11 0.34 



2 2 

1.674465 
0.035 



3.599230 
0.19 



-5 -> +| 
4.565917 
0.17 



+1 ~> -5 
6.387734 

0.155 



As a good approximation, calculating the dressed en- 
ergies through the effective Lande g-factor of Eq. (J2|) we 
get the following expression of the magic B m field as a 
function of the uj r f and x m parameters: 



4 2) 



m F'9pS{^x m ) - m F g s S{x m ) 



(8) 



Table I reports at fixed ui r t few magic pairs (x m , B m ), 
at increasing B m values, for n polarization atomic clock 
of the fcrmionic species of present interest. The uj r f /2tt = 
2 kHz choice, producing a consistent table for all the 
fcrmionic atoms, leads to very small magic field values for 
the mercury atom. Higher values, more easily manage- 
able in the laboratory, are simply obtained by increasing 
the rf frequency and applying the Eq. (JSJ scaling. 

The magic field values are calculated at a fixed Cl/u> r f 
value, that implies a very large precision in the setting of 
the fl and u> r f parameters. For a practical application, 
the stability of those quantities becomes an important 
issue. While a very high stability of the rf frequency is 
not a problem, the f2 accuracy could be an issue. We 
have explored the Q stability required in the operation of 
an optical clock, and the results are reported in Fig. @] 
A change in 57, or precisely a change in x by one part 
in ten thousand corresponds to a fractional shift of the 
optical clock at the 10" 17 level. In order to reach the 
ultimate limit of the alkali-earths optical clocks be- 
cause only the a; ratio of the rf quantities is important for 
the rf engineering, the ft variations may be compensated 
by acting on the rf frequency. Thus, a feedback on the 
rf frequency should reach the required ratio stability, in 
order to reach a 10~ 20 level in Fig. 01 In practice the rf 
stability could be matched to the actual accuracy of the 
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FIG. 4: (Color online) Fractional clock shift, measured in 
10 -20 relative units, for the 87 Sr rriF = —9/2 it transition, 
for a rf amplitude fluctuation by 0.01 percent (left scale) and 
for correlated variations of the same amount applied to both 
amplitude and rf frequency (right scale). Operating point 
x r f = 0.96... and UV//2-7T = 2 kHz, leading to the magic 
B — 2.06 fiT magic value. 



optical clock. For large uj r f excursions, the efficiency of 
this compensation is limited by the w 2 ^ dependence in 
the B m numerator of Eq. ([5]) . 

We have verified that the shift produced by virtual 
transitions induced by the rf field between the 3 P fine 
structure levels is negligible compared to the finally 
aimed clock accuracy. The presence of vectorial ac shift 
contribution to the atomic energies introduces an ad- 
ditional shift of the clock levels. When the ac shift is 
comparable to the Zecman shift, the key features of the 
rf dressing, crossing nodes, compensations of the clock 
shift and magic values, are obtained also for this case, 
with "magic" (x m ,B m ) pairs depending on the specific 
ac shift. Schemes for the rf compensation of the ac tenso- 
rial part appearing in the 87 Sr optical clock only should 
be investigated. 

The combination of well-engineered optical-trapping 
potential and of rf quantum engineering represents an im- 
portant tool in the investigation of alkali-earth clock sys- 
tematics. We have verified our scheme feasibility within 
the operation regime of the present optical clocks. The 
averaging over the Zeeman components of the optical 
clock transition is directly performed by the rf dressing 
of the atomic system. For an implementation within an 
optical lattice where magnetic fields are created synthet- 
ically [15| . the dressing magnetic field may be originated 
by the rf modulation of the lattice depth, at least for fre- 
quencies low enough for an atomic adiabatic following. 

The synthetic rf controlled magnetism may be applied 
to other atomic/molecular and solid-state physics con- 
figurations. Beside compensating the residual Zeeman 
contribution to a super stable optical clock based on a 
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nuclear transition p9j , the rf engineering may be applied 
to design an artificial quantum transition with specific 
Zeeman properties, as a two-level superconducting sys- 
tem driven by an oscillatory field [30| , and to the control 
of spin coherent dynamics and transport in semiconduc- 
tor systems [3l|, [HJ ■ 
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